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Abstract 

Q H ' We construct the general solution of the Yang-Baxter equation which is represented as an integral 

operator with an elliptic hypergeometric kernel acting in the space of functions of two complex variables. 
It intertwines the product of two standard L-operators associated with the Sklyanin algebra, an elliptic 
deformation of s£(2). This R-operator is constructed from three basic operators Si, S2 and S3 generating 
the permutation group of four parameters 64. Validity of the key Coxeter relations is based on the elliptic 
beta integral evaluation formula and the Bailey lemma associated with an elliptic Fourier transformation. 
' The operators Sj are determined uniquely with the help of an elliptic modular double. 
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1 Introduction 



The Yang-Baxter equation (YBE) 

R12O - v) R 13 (u) R 23 (v) = R 23 (u) Ri3(«) Mi 2 (u - «) (1.1) 

plays a key role in the theory of completely integrable quantum systems [IHZ]- Its general solution 
is described by the operator Rj^(u) acting in the tensor product of three (in general different) spaces 
Vi <S> V2 <8> V3 . The indices i and k show that Rjfc (it) acts nontrivially in the tensor product Vj <S> 
and it is the unity operator in the remaining part of Vi (8> V2 ® V3 . The operator Rjfc(u) depends on 
the complex spectral parameter u and is called the R-matrix. 

For spaces Vj of finite dimension there are three increasing levels of complexity of known YBE 
solutions described by rational, trigonometric and elliptic functions. Investigation of the most compli- 
cated elliptic level was initiated by Baxter [8] for the case when all Vj-spaces are two-dimensional. In 
a more general setting, when one of the spaces becomes infinite-dimensional, a major role is played 
by the Sklyanin algebra [9,:1Q]. Our main goal consists in the construction of general solution of the 
Yang-Baxter equation when all spaces Vj are infinite-dimensional and Rjfc(it) are described by integral 
operators. In this case the hierarchy of solutions of the YBE is attached to the plain hyper geometric, 
hypergeometric and elliptic hypergeometric functions in the increasing order of complexity. An 
explicit realization of the R-matrix as an integral operator in the simplest situation of symmetry algebra 
s£(2) has been considered in detail in [12]. In the present work we discuss only R-operators related to 
the most complicated elliptic level. 

Elliptic hypergeometric integrals were introduced in |13pi4] . They define the general form of elliptic 
hypergeometric functions which cannot be approached by infinite series [15] because of the convergence 
problems. The discovery of these functions and various relations for them formed a breakthrough 
in the theory of special functions. These functions generalize all previously constructed functions of 
hypergeometric type and they still obey the properties characteristic to classical special functions [16]. 
In particular, elliptic beta integrals |lipi3] form a new class of exactly computable integrals. A kind of 
elliptic Fourier transform was introduced in |17] as an integral generalization of the well known Bailey 
transformation |16] . An elliptic extension of Faddeev's modular double [18] was introduced in [19] . Both 
these constructions will play a major role in our consideration below. 

We use the general strategy of building integral operator solutions of YBE whose initial steps were 
discovered in [20J. Its formulation was completed at the rational level in [21|,l22j . where an SL(N,C)- 
invariant solution of YBE related to ^4 n -root system has been constructed. In [23J this method was 
used also at the elliptic level employing some formal infinite series. Here we apply it for building 
integral operator solutions related to elliptic hypergeometric integrals. First we define some useful formal 
operators Si , S2 and S3 performing permutation of parameters in the defining RLL-relation. One of these 
operators is an intertwining operator of the Sklyanin algebra. Then we build these operators explicitly 
as integral operators with elliptic hypergeometric kernels. Finally we prove the Coxeter relations for 
these operators, for which the elliptic beta integral and the related integral Bailey lemma play a crucial 
role, and confirm that they indeed generate the 64-group. As discussed in the concluding section, our 
results have applications to an interplay between integrable 2d spin systems and 4d super symmetric 
gauge field theories. 
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2 Sklyanin algebra 

In the simplest case of equation (jl.ip all the spaces are two-dimensional, Vj = C 2 . Then the most 
general R-matrix has the form 



1 / 6a+l{u + 77 ) 
-12 = 77/, W a{ u ) a a ® °"a 5 w a{u) = — p— , (2.1) 



2^~»— " ' -v-/ ^ +i(T?) 

where <r a are the standard Pauli sigma-matrices and we use the shorthand notations 9 a (x) = 9 a {x\T) for 
the standard Jacobi theta-functions with modular parameter r. The definitions and useful formulae for 
# a -functions are listed in the Appendix. This R-matrix was found by Baxter in solving the eight-vertex 
model [21E] 

At the next level of complexity of relation (jl.ip one of the spaces, say V3, is arbitrary, and Vi, V2 = 
C 2 . In this case the R-matrix R^-u) (and R23(u)) is known as the quantum L-operator or the Lax 
matrix. It is a matrix with two rows and two columns 

/ \ w •> 1 f \ „„ 1 ( wq(u) S° + ws(u) S 3 wi(u) S 1 — iwo(u) S 2 \ , n „. 

L13 (u) = L(«) := -\w a (u)a a ® S a = - uv ' x . , \ 2 1V \ 2V { 3 , 2.2 
2^ 2 \ wi(u) b + iW2{u) S Wo(u) S — W3(u) a J 

where the matrix element entries S a are some operators acting in V3. In this case the equation for 
L-operator is the Yang-Baxter equation of the form 

R 12 (u - v) L n (u) L 23 (v) = L 23 (v) L 13 (u) R 12 {u - v) , (2.3) 

where R12 (it) is Baxter's R-matrix (|2. 1 j) . This equation is equivalent to the following set of relations for 
four operators S , S 1 , S 2 , S 3 forming the Sklyanin algebra [9lll0j: 



S° S° - S a S° = i J^ 7 • (V S 7 + S 7 S' 3 

where the triplet (a,/3,7) is an arbitrary cyclic permutation of (1,2,3) and the structure constants J^ 7 
are parameterized in terms of theta functions as 

One can write J a p = J/3 j J ° , where 

_ g 2 (2ry)g 2 (0) , _ e 3 (2 V )9 3 (0) _ e 4 (2 V )9 4 {0) 

There are two Casimir operators commuting with all generators: [Ko, S a ] = [K 2 , S a ] = 0, 

3 3 



K = S a S a ; K 2 = ]T J a S a S a . 



a=0 a=0 

We shall use the explicit realization of generators as difference operators [9j[T0] 

[ga $] (z) = WH+iML a+i (2z _ . $(z + v) _ 9a+i { _ 2z _ 2rj£) . Hz _ T}) , 2 - , 
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In this realization the Casimir operators reduce to the following scalar expressions 

Ko = 40?((2* + l)»j) ; K 2 = 40i(2(* + l)»7)0i(2£»7). 

The variable £ is called the spin. It labels the algebra representations since it fixes (together with r\ and 
r) the Casimir operator values. 

There exists a useful factorized representation for the L-operator 

j, , 1 ( 3 (z- Ul ) -d 3 (z + Ul ) \ ( \(h(z + u 2 ) 9 3 (z + u 2 )\ 

MUl,U2j e 1 (2z)\-9 4 (z-u 1 ) h{z + u x ) ) \ e-* d * )\h{z-u 2 ) 9 3 (z-u 2 )J> 

where e^ 2 is a shift operator, e 719 * f(z) = f(z+rj). New parameters u\ = ^+r](£+^) and u 2 = ^—rj(£+^) 
are simple linear combinations of the spectral parameter u and the parameter £ labeling representations 
of the Sklyanin algebra. Here we use notation 9 a (x) = 9 a (x\^) for theta- functions with the modular 
parameter J. 

In the case of integer or half-integer £ there exists (21 + l)-dimensional subspace 0^ of even theta- 
functions of order A£ which is invariant under the action of generators S a . 

The L-operator (|2.2p is not unique. For instance, the operator <7;,L, where 07, is any Pauli matrix, 
is also a solution of equation (12. 3j) . This follows from the fact that Xf, := Ob <S> (Jb obeys the properties 
X? = 1 and XbWLi>z(u)Xb = Ri2(it). This freedom leads also to the existence of nontrivial automorphisms 
of the Sklyanin algebra |10| . 

The top level of complexity of the R-matrix corresponds to the most general situation when all spaces 
Vj are infinite-dimensional. In this case one deals with the general solution of Yang-Baxter equation. 

In the next section we explain step-by-step our general strategy for building this solution, which 
is essentially the same as in [23] where the important role is played by the intertwining operator. For 
integer or half-integer £ such an intertwining operator was constructed first by Zabrodin [24] as a finite 
sum of the powers of the finite-difference operator e~^® z . Its straightforward extension to arbitrary 
values of £ as an infinite series proposed in [25J has only a formal meaning due to the convergency 
problems. Nevertheless, the needed Coxeter relations were verified in [23] by checking the equality of 
coefficients in two such formal infinite series with the help of the Frenkel-Turaev summation formula [15] . 

Here we put the construction on a firm mathematical ground by using a different general ansatz for 
the intertwining operator which is more useful for practical applications. The key ingredients needed 
for the completion of this program is the elliptic beta integral [111113] . the most general known exact 
integration formula generalizing the Euler beta integral, and the elliptic Fourier transformation of [T7] 
which is defined precisely with the help of the needed intertwining operator. 



3 General construction 

We solve YBE (jl.lh when all spaces V» are infinite-dimensional in two steps: 

• on the first stage, we solve a defining RLL-relation using as elementary building blocks some 
simple operators Si, S2 and S3. The key structural entries at this step are Coxeter relations for 
the operators Si, S2 and S3 validity of which is guaranteed by the elliptic beta integral evaluation 
formula; 

• on the second stage, we prove that the operator Ri2(u) found from the defining RLL-relation 
obeys the general YBE. 

Consider a realization of the YBE different from the previous ones, namely, when the spaces Vi and 
V2 are arbitrary and the space V3 is two-dimensional. Then equation (jl.lh is reduced to the defining 
equation for general R-matrix called RLL-relation [9|I26]: 

Ri 2 (u - v) U(u) L 2 (v) = L 2 (v) U(u) R 12 (u - v) . (3.6) 
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Here we use compact notation: the index k in indicates that the generators S a entering this matrix 
are the operators acting in the space Vfe, i.e. : Vfc — > Vfe. The Sklyanin algebra generators Sf 
and S2 act in different spaces and evidently commute, [Sf,S2] = 0. Matrices L/% in equation f|3.6[) are 
multiplied as usual 2x2 matrices. 

Due to the non-uniqueness of representation of the L-operator, there are several possible forms of 
equation (|3.6[) with different R-operators labeled by two indices a and b enumerating possible L-operators 

M.f 2 (u - v) a a Li(«) a b L 2 (v) = a b L 2 (v) a a U(u) M.f 2 (u - v) . 

For a technical reason, which will be clear a posteriori, we fix a = b = 3 from the very beginning. In 
this case it is possible to cancel one 03 and our main defining RLL-relation takes the form 



KL2 



(u - v) Li(it) 03 L 2 (v) = L 2 (v) 03 Li(n) R 12 (u - v) . (3.7) 



We assume that Vi is the space of functions of the complex variable z\ and V2 is the space of functions 
of the complex variable z 2 . Respectively, the space Vi ® V2, where M12 is acting, is the space of functions 
of two independent variables z\ and z 2 . 

It is convenient to extract from the R-matrix the permutation operator Ri2(u) := Pi2Ri2( M )> where 
the permutation operator interchanges arguments P12 ${zx, z 2 ) = &(z 2 ,zi). The defining equation for 
the operator R12 has the following explicit form 

Ri 2 (u - v) Li(m,n 2 ) o- 3 L 2 (vi,v 2 ) = Li(vi, v 2 ) cr 3 L 2 (u 1 , u 2 ) Ri 2 (u - v) , (3.8) 

where the operators z,d z in the Sklyanin algebra generators (|2.5p entering Li are replaced by z\,d zi 
and in L2 they are replaced by z 2 , d Z2 . We use also the following notation 

«x = (4 + 5) >"* = %-v(*i + l) ; * = h*fa + l)> va = h''( ia + ii) (3 - 9) 

for combinations of spectral parameters and spin variables. For further needs it is convenient to assume 
that these parameters do not depend rj and r (i.e., to assume that spectral parameters u, v and g\ j2 = 
7/(2^2 + 1) are independent on r\ and r). 

Equation (13. T[) admits a natural interpretation: the operator R12 interchanges the set of parameters 
(ui,u 2 ) from the first L-operator with the set of parameters (^1,^2) in the second L-operator. It is 
useful to combine these four parameters in one set in the natural order u = (ui, u 2 , v%, v 2 ). Then the 
operator 

R12C" -v) = R i2 (u) = Ri 2 (u 1 ,u 2 \vi,v 2 ) 
corresponds to a particular permutation s in the group of permutations ©4 of four parameters: 

s -)■ Ri 2 (u) ; su = s(ui,u 2 ,vi,v 2 ) = (vi,v 2 ,u 1 ,u 2 ). 

Any permutation from the group 64 can be composed from the elementary transpositions si, S2 and S3: 

SlU = (u 2 ,Ui, Vi,V 2 ) , s 2 u = (ui,v 1 ,ui,v 2 ) , S 3 U = (til, U2,V 2 ,Vt), 

which interchange only two nearest neighboring elements in the set (ui,u 2 ,vi,v 2 ). For example, the 
permutation s has the following decomposition s = s 2 siS3.s 2 . It is natural to search for the operators 
Si(«i, u 2 , vi, v 2 ) = Sj(u) representing these elementary transpositions in L-operators 

Si S3 S2 

(uT^m,vT^2) ; (Ul ,U2?Ol, v 2 ). 
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These operators obey the following defining relations 

Si(u)Li(ni,u 2 ) = Li(u 2 ,«i)Si(u) ; S 3 (u) L 2 (v 1 , v 2 ) = L 2 (u 2 , vi) S 3 (u), (3.10) 
S 2 (u)L 1 (ui,W2)<T3L 2 («i,U2) = Li(ui,vi)ct 3 L 2 (-u 2 ,v 2 )S 2 (u) . (3.11) 

Our first step consists in direct construction of these operators (see the next section). Having these 
operators we can construct the R-matrix. 

Proposition 1. The operator Ri 2 (u), 

R i2 (u) = S 2 (sis 3 s 2 u) Si(s 3 s 2 u) S 3 (s 2 u) S 2 (u) (3.12) 
satisfies equation \3.8\) as a consequence of relations \3.10\) - L3.11\) for the operators Sj. 

The proof reduces to a direct check, which is quite simple. This expression for the R-matrix corre- 
sponds to a special decomposition of the permutation s: s = s 2 sxs 3 s 2 . We will see that operators Si 
depend on their parameters in a special way 

Si(u) = Si(«i - u 2 ) ; S 2 (u) = S 2 (u 2 - vi) ; S 3 (u) = S 3 (vi - v 2 ), (3.13) 

so that the operator Ri 2 (u) depends on the difference of spectral parameters u — v as it should 

R\ 2 {ui,u 2 \vi,v 2 ) = S 2 («i - v 2 ) Si(ui - vi) S 3 (-u 2 - v 2 ) S 2 (u 2 - vi) . (3.14) 

We have thus the correspondence 

Si — ► Sj(u) ; SiSj — > Sj(sju) Sj(u). (3.15) 

In order to prove that we have a representation of the permutation group ©4 it remains to prove the 
defining Coxeter relations for the generators 

SiSi = 1 — > Sj(sju) Si(u) = 11 ; sis 3 = s 3 si — > S x (s 3 u) S 3 (u) = S 3 (siu) Si(u) , (3.16) 
sis 2 si = s 2 s\s 2 — ► Si(s 2 siu) S 2 (siu) Si(u) = S 2 (sis 2 u) Si(s 2 u) S 2 (u) , (3.17) 
s 2 s 3 s 2 = s 3 s 2 s 3 — > S 2 (s 3 s 2 u) S 3 (s 2 u) S 2 (u) = S 3 (s 2 s 3 u) S 2 (s 3 u) S 3 (u) . (3.18) 

The proof of all these relations will be given in Sect. [5] and the key tool needed for that is the elliptic 
beta-integral [13J . 

Proposition 2. The operator Pi 2 Ri 2 (u) satisfies the Yang-Baxter equation (jl.ip as a consequence of 
the Coxeter relations f)3. 16|) - (|3. 18|) . 

Consider the group of permutations of six parameters ©6 i n the product of three L-operators. The 
following permutation of L-operators 

Li(ui,u 2 ) a 3 L 2 (v 1 ,v 2 ) a 3 L 3 (wuW2) -> Li(wi,w 2 ) cr 3 L 2 («i, v 2 ) a 3 L 3 («i, u 2 ) 

can be realized in two different ways as shown on the figure below 
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R23{ui,U2\wi,W 2 ) 

L 1 (vi,V2)a 3 L 2 (u 1 ,u 2 )cr 3 L 3 (wi,W2) Li(«i, v 2 )cr 3 L2(wi, w 2 )a 3 L 3 u 2 ) 



Ri 2 {u 1 ,u 2 \v 1 ,v 2 ) 



L 1 (m,u 2 )a 3 L 2 (vi,v 2 )a 3 L 3 (w 1 ,w 2 ) 



R23(vi,V2\wi,W 2 ) 



L 1 (ui,U 2 )cT3L2{w 1 ,W2)a 3 L 3 (vi,V2) 



Rl2{vi,V 2 \wi,W 2 ) 



Li(i«i , w 2 )a 3 L 2 (vi, v 2 )cr 3 L 3 (u 1 , u 2 ) 



R23(UI,U 2 \VI,V 2 ) 



R 12 (ui,u 2 \wi,w 2 ) 



Li(i«i , w 2 )a 3 L 2 (ui, u 2 )cr 3 L 3 (v 1 , v 2 ) 



The Yang-Baxter relation for the operator Ri2(u) = Ri 2 (ui, u 2 \v\, V2) coincides with the condition 
of commutativity of this diagram which yields 

Ri2{vi,v 2 \w 1 ,w 2 )R 23 {ux,U2\w 1 ,w 2 )R 12 (u 1 ,u 2 \vx,v 2 ) 

= R 2 z(u 1 ;u 2 \v 1 ,v 2 )R 12 (ui,u 2 \w 1 , w 2 )R 23 (v 1 ,v 2 \w 1 ,w 2 ) . (3.19) 
Introducing the needed permutation operators Fjk into this relation one comes to the original equation 

4 Elementary transpositions and intertwining operators 

We shall use the factorized form of the L-operator which allows one to simplify considerably all calcu- 
lations 

1 / e^ 9 \ 

L(u 1 ,u 2 ) = Q^z) ' M<yZ ~ " i; Z + ' i e _r?9 J • N ( 2:_u 2;^ + u 2 ), 

where 

This factorization of the Lax matrix is based on the addition formula 

#4 (a) h (b) - h (b) 6 3 (a) = 2 0i(a + 6) 0i(a - 6). (4.2) 
The product of matrices ./V and M has the form 

N(« i;i)1 ).Mfe fe ) = 2.f «l(h-«2)h(h+«2) -h(h-WAh+h) \ : (4 . 3) 

V 6»! (ai - a 2 ) 6»i (ai + a 2 ) -6>i (ai - f> 2 ) 6>i (ai + b 2 ) J 

in particular, 

N{a;b)-M(a;b) = -2-6 1 {a-b)9 l (b + a) ( J M. (4.4) 

To avoid lengthy formulae we use compact notations N(a=F&) = N(a — 6; a+6), M(a=Ffr) = M(a — b; a + b) 
and (a =F 6) = 0\ (a — b) 8\ (a + b) , so that 

1 / e^ \ 
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and 



N(o i; b x ) ■ M(a 2 ;b 2 ) = 2\e l{aiT a2) _ ex \ ai Tb2) )~ (4-5) 
Consider first the denning relation for operator S2 (|3.1ip . 

/ e vdi \ / e r > d2 \ 

S 2 M(ziT«l) ( e - v di J N{z l ^u 2 )a 3 M(z 2 ^v 1 ) I Q ^ J N(z 2 =F v 2 ) = 

( e'l 91 \ / e 1 ! 92 \ 

= M(zi=f«i) ( e - v a 1 J N(«i T vi) 0-3 M(z 2 T 1*2) ( Q j N(z 2 t w 2 ) S 2 . 

We underlined the matrices which can be canceled under the commutativity condition [S 2 , z\\ = [S 2 , z 2 \ = 
0. This observation suggests that S 2 is just a multiplication operator: 

[S 2 $] (zi,z 2 ) = S(z u z 2 ) ■ $(zi, z 2 ) . 

Consequently, the operator S 2 commutes with the matrices M(zi =F u\) and N(z 2 =F v 2 ), so that they 
both cancel from the equation and we obtain a much simpler defining relation for the function S(zi, z 2 ): 

S(zi,z 2 ) I Q e _ r?9l J •N(ziTn 2 )o- 3 M(z 2 =Fwi) I Q g _^ 2 J = 

/ e »?9i q \ / e rid2 q \ 

= 1 Q J •N(2iTt)i)(73M(22TM2) I Q ^ v g 2 J S {z X , Z 2 ) , 

or in the equivalent form 

S(zi - 77, z 2 ) \ , T , \ i\ t / \ 

n ar , \ N zi =F« 2 )cr 3 M z 2 =F = 

SOi +r),z 2 ) J y 



= N(z 1 ^vi)a 3 M(z 2 ^u 2 ) 

Note that 



S(zi,z 2 +77) 

S(zi,z 2 -r?) 



/ e 27ri( 2l + 5+ M2 ) \ 

N(^=F« 2 )o-3 = -^ Q e2li(2li _ tt2) jN( Zl + i T « 2 ), 

and then using the general formula for the product N(ai, b x ) ■ M(a 2 , b 2 ) we find 

/ p 27ri( 2 i + f +u 2 ) n \ 

N(z x =F u 2 ) a 3 M(z 2 T v x ) = -2 Q e2m(2l+f _ U2 ) J ■ 

/ 6»i (zi + u 2 + § =F O2 - vi)) -9i (zi + u 2 + § =F (z 2 + vi)) \ 
' \ 6»i (zi - u 2 + § =f (-22 - ui)) 6*i (z 2 + v x + § =f (21 - u 2 )) J ' 

The second matrix N(zi =F i>i) (73 M(z 2 =p u 2 ) is obtained after the simple permutation u 2 <^ v\ 



N(zi =F vi)o- 3 M(z 2 =fu 2 ) 



e 27ri(zi + 5+i)i) Q 

q e 27ri(zi+5-m) 



?1 (zi + + | =F (z 2 - u 2 )) -9 X (zi + «i + \ =F (z 2 + n 2 )) 
?i (zi - + § =F (z 2 - u 2 )) 81 (z 2 + u 2 + \ =F (zi - i>i)) 

Note that the terms depending on u 2 + v\ cancel in the final equations. After all these manipulations 
we come to the following system of four finite-difference equations 

e 2 7 ri( U2 -« 1 ) q^ Zi + Z2 + U2 _ ^ + r) g( Zl _ ??) Z2 ) = +Z2+Vl _ U2 + r- ) S ( Z1) Z2 + r? ) ; 
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e - 27 ri( W2 - Vl ) +Z2 _ U2 + Vl + r^ g (zi + r?j = + ^ _ Wi + ^ + r j g^ % _ ^ > 

e 2 " 1 ^"" 1 ) - z 2 + u 2 - Wl + §) S(z! - 7?, z 2 ) = 9 1 {z 1 - z 2 + v t - u 2 + §) S(*i, z 2 - r?) , 

e -2*i(t«-«i) ^ _ ^ _ ^ + ^ + S(zi + ?7j Z2 ) = 0^ _ Z2 _ Wi + U2 + r j g^ ^ + ^ _ 
Their structure suggests to look for the solution in the factorized form 

S(z 1 ,z 2 ) = <5> + (z 1 + z 2 )-$^(z 1 - z 2 ). 
Then, in each equation one of the ^i-factors drops out and we obtain equations 

e 2ni{u2 - Vl) 0i (z + u 2 - vi + |) $±(z - 7?) = 0i - u 2 + ui + §) $±(* + ry), (4.6) 
or in the equivalent form 

* ± (* + 2,) = e-(—) fli(* + t*2-*i+q+5) . $ ±{z y (4 . 7) 

A particular solution of these equations is described by a ratio of elliptic gamma functions T(z) = 
r(z|r, 2rf) (see the Appendix), 

+ _^>*£M.. Wi . W -^M. (4,) 

Note that this solution $(z) is well defined only for Im(r/) > 0, whereas the initial equation does not 
assume such a restriction on the parameter rj. A solution valid for Im(ry) < has the form 

= r(z + b-2rj\T,-2rj) 
1 ' T(z + a-2r]\T,-2r))' 

Solutions valid for Im(ry) = are described in terms of the modified elliptic gamma function |1HI14| 
which we do not describe here. For simplicity, we stick to the region Im(^) > 0. 
Using this result we can write the general solution S(zi,z 2 ) in the form 

T{z 1 + z 2 + u 2 -vi + 7] + ^)T(z 1 -z 2 + u 2 -v 1 +7 ] + ^) 

S{Zl,Z 2 ) = — ■ ■ ■ — F -— ■ ■ — pr^l,^), 4.9) 

T(Zi + z 2 - U 2 + Vl + T) + \) Y\Z\ - z 2 - U 2 + Vl + 7] + \) 

where tp is an arbitrary function satisfying the periodicity conditions 

<p(zi + 2r],z 2 ) = <p(zi,z 2 + 2rf) = <p(zi + rj,z 2 + rj) = (p(zi,z 2 ). 
Using reflection formula for the elliptic gamma function, T(x) T(—x + 2rj + r) = 1, and the notation 

T (=fx T y + a) = T (—x — y + a) T (—x + y + a) T (x — y + a) T (x + y + a) 
it is possible to rewrite the latter formula in a much more compact form 

S(zi, z 2 ) = T (tzi Tz 2 + u 2 -v 1 +r] + ^) (p(zi,z 2 ). 

The functional freedom <p(z\,z 2 ) strongly influencing the final results will be fixed in the section 
dedicated to the elliptic modular double. For the moment we simply set ip = 1. We remark also that 
our choice of the L-operator in the form 03 L is done for technical reasons in order to avoid unpleasant 
exponential factors in the expression for S 2 -operator emerging for other choices. 

Let us consider now defining equations (|3.10p for operators Si and S3. Permutation of the parameters 
u\ = I + 77 (^1 + and u 2 = ¥■ — rj (£\ + ^) is equivalent to the change of the spin l\ — > — 1 — 1\ and 
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similarly the permutation of parameters v\ = | + rj (£2 + \) and V2 = § — rj (£2 + 5) is equivalent to the 
change £2 — > — 1 — ^2- In the Lax matrix ()2.3[) only generators S a depend on the spin, therefore defining 
equations A3. 1QH can be rewritten in terms of the S a -generators only: 

Si • S^z^h) = S a ( Zl , -1 - h) ■ Si , S 3 • S a (z 2 ,£ 2 ) = S a (z 2 , -1 - £2) ■ S 3 • 

The meaning of these relations is the following: the operator Si intertwines representations with the 
spins £\ and —\ — £\ realized in the space of functions of the variable z\ and the operator S3 intertwines 
representations with the spins £2 and — 1 — £2 realized in the space of functions of the variable Z2- Note 
that in terms of the variable g = r/(2£ + 1) this corresponds to the simple sign change g —> —g. 

Evidently the operators Si and S3 are equivalent to each other differing only by the spaces of 
representation. Let us construct the general intertwining operator W defined in two equivalent ways: 
either as a solution of the matrix equation for the L-operator 

W • L(ui,u 2 ) = L(u 2 ,ui) • W, 

or, alternatively, as a solution of the system of equations involving generators of the Sklyanin algebra 

W- S a (z,£) = S a (z,-1 -£) ■ W. (4.10) 

For integer or half-integer £ such an intertwining operator W was constructed first in [24] as a finite 
sum of the powers of the finite-difference operator e~ rtdz (see below). A formal extension of this W to 
infinite series and arbitrary values of £ proposed in |25j does not represent a well defined operator. 

In this paper we use a different approach inspired by the elliptic hypergeometric integrals discovered 
in [13|J14|. Namely, we construct the intertwining operator as an integral operator 

[W I dxn(x)W(z,x)$(x) ; (4.11) 

J 

= r(2 C x)^-2x) ' c = 5 (e4 "" ; e4 ""' ) ~ C^i^- 

The measure fJ,(x) is chosen in a special way as a solution of the relations 

Vi(2x) 0i(2x) 

and the choice of normalization C is dictated by the elliptic beta integral [13] as will be evident below. 
With this choice of the measure the generators of Sklyanin algebra have simple conjugation rules 

f dxn(x)V(x) ■ [S a (x,£)<S>] (x) = (-l) Sa '° [ dxfj,(x) [S a (x,-1 (x) ■ $(x), (4.12) 

Jo Jo 

or in the compact form 

[S a (^)] f = (-) Sa >° S a (-l-£). 
The L-operator has the form with the weights W{(u) as described above in (|2.2|) . 

L(«) = y> a («)o«®S-=( ^St^lt ^"^IfaV (4-13) 
K J ^ \ W!(u) S 1 + iw 2 (u) S 2 w (u) S° - w 3 (u) S 3 J v ; 

Assume for the moment that w a (u) = 9 a +\{u)/9(r]), which differs from our initial choice by the shift 
u — y u — fj. Then the function wq(u) is odd: wq{— u) = —wq{u) and functions w^{u) are even for 
k = 1, 2, 3 : Wk(—u) = Wk(u). Therefore the conjugation rule for the whole L-operator is 

[L(u,£)]i = L(-u,-l-e), 
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or in another parametrization [L(iti,tt2)] = L(u 2 ,ui), where the parameters u± and 112 differ from u\ 
and U2 by the sign of the spectral parameter: 

ui = -+r}£; u 2 = - -r){l+ I) ,u\ = -- + 7]£ ; u 2 = -- - 77 (£ + 1) . 

After these preliminary considerations we rewrite relation (|4.10p as an equation for the kernel of 
operator (|4.1ip 

L(u 2 , ui) W(z, x) = L(u 2 ,ui)W(z, x). (4.14) 

To avoid misunderstanding note that all matrix elements of L(u2,mi) and L(u2,«i) are finite-difference 
operators acting on the variables x and z, respectively. Using the factorized form of the L-operators we 
obtain 



M(x T « 2 )( Q e . vdx J N(xtui)W(z,x) 



1 _ s ( e vdx 

01 (2x) 

1 / e V&z Q 

M{ztu 2 )[ n I N(zT!ii)W(z,4 



This equation is similar to the equation for function S 2 (zi, z 2 ) from the previous section, so we proceed 
as before. First, we multiply by the matrix N(x =F U2) from the left-hand side 

-^—^ • N(x T Sa) M(x T «2) ( e \ ) N(x T «i) W(z, x) = 

= • N(x T «2) M(z T «a) ( ^ e \ ) N(* =F ui) W(z, x) 

and use auxiliary formulae (|4.3p and (|4.4p 



^i(2^)-( e ^ e \ )-N(x T n 1 )W(^,x) = ^ 
01 (x + u 2 =F (2 - u 2 )) —6\ (x + u 2 =F (z + 1*2)) 



?i (x - u 2 T (z - us)) -9i(x-u 2 T(z + u 2 )) J\ e"^ 



•N(z=fui)W(z,x). 



Next we multiply by the matrix M (a ^ (77 + ui)) from the right-hand side and use the relation which 
can be easily checked 

C e -*»« ) N(Z T Ul) M ( ° T {V + Ul))la ^ = 201 {2Ul + 2V) dl {2Z) ( & 7 e"°^ ) • 

The notation M (a =F (v + ' li i))| a _ >2 means that we substitute a — > z after acting by the operators e^ 2 . 
More precisely, we multiply from the right by the matrix M (a =p (r/ + u\)) independent on z, compute 
matrix products and move all operators e^ 2 to the right and only after all these manipulations we 
set a = z (which is allowed since a is an arbitrary parameter). This procedure is equivalent to taking 
initial equations and forming from them some linear combinations using as coefficients matrix elements 
of M(*=F(»7 + ui)). 

Because 9\ {2u\ + 2rj) = — 0i(2u 2 ), one obtains the product N(x =F u\)M (z =p (7/ + ui)): 



e vdx \ / 0i (x + ut ^ (z - T) - ui)) -0i (x + ni T (z + rj + ui)) 
e~ vdx ) \ 0i (x - u\ T (z — r] — u%)) -0% (x-ui^(z + r] + ui)) 

-0i (x + u 2 =F (z - u 2 )) 0i (x + n 2 =f (z + u 2 )) \ / , 

-0i (x - u 2 T {z - u 2 )) 0i (x - u 2 T {z + u 2 )) )'\ ' ' 



W 
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There are terms depending on u\ — u\ = U2 — U2 = u. However, the operator W does not depend on 
the spectral parameter and one can expect that such terms cancel, and, indeed, this cancellation can be 
easily checked. Using relations u\ + u\ = 2r\i , U2 + U2 = —2r](£ + 1) and u\ — U2 = f](2t + 1), we arrive 
to four finite-difference equations 

-6i(x + z - ui + u 2 + rj) W(z, x - rj) = 9\{x + z + u\ - u 2 + rj) W(z + r/,x) , 

—9\{x + z + u\ — U2 + rj) W(z, x + rj) = 8\(x + z — u\ + 112 — rj) W(z — 77, x) , 
—9\{x — z + u\ — U2 + rj) W(z, x + rj) = 6\{x — z — u\ + U2 — rj) W(z + r],x) , 
—9\{x — z — u\ + U2 — rj) W(z, x — rj) = 9\{x — z + u\ — U2 + r]) W(z — rj, x). 
Let us search for solutions in the factorized form 

W(z, x) = + z) ■ - z). 

Then in each equation one of the $±-factors drops out and we come to difference equations 

-e 1 (z-u 1 + u 2 - rj) <£>±(z - rj) = 9i{z + m - u 2 + rj) $±(z + rj), (4.15) 

or in the equivalent form 

* ± (* + H = ~ ^ Z - Ul + U "\ , • * ± W . (4.16) 

9 l {z + u l -u 2 + 2r]) ±K ' y 1 

Assuming again that Im(r/) > 0, one finds a particular solution 

_iri 2 

*±{z) = — \ \ , a = Ul -u 2 . 4.17 

T(z + a + 2t] + t) 

Finally, the general form of the intertwining operator kernel is 

, , e " v T(x + z-ui+u 2 ) T(x- Z-U1+U2) , . , , 

W(Z,X) = — 7— r— -(p(Z,x), (4.18) 

v ' ; T{2u2-2u l )T(x + z + u l -U2 + 2r ] + T)T(x-z + u l -U2 + 2r, + Ty y ' " v ' 
or in compact notations 

_S.( 2 2 +a; 2) 

W(z, x) = — — _ — T(^x T z - ui + u 2 )ip(z, x), (4.19) 

where (p(z + 2rj, x) = <p(z, x + 2rj) = ip{z + rj, x + rj) = ip(z, x) is an arbitrary periodic function. The 
latter functional freedom will be fixed in the section dedicated to the elliptic modular double. Similarly 
to the previous section, here we simply set ip = 1. Note that before fixing this normalization we singled 
out an elliptic gamma function multiplier T{2u2 — 2u\) in the denominator. 

This result was obtained for u — )■ u — rj shifted form of the L-operator. To obtain the intertwining 
operator for the original L-operator (|2.2j) we shift u — > u + rj in (14. 191) . which is simply equivalent to 
returning to the notation (|3.9p for parameters. Actually the intertwining operator does not depend on 
the spectral parameter at all since it involves only the difference u% — U2- 
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5 Coxeter relations and the elliptic beta integral 

In this section we prove that the derived operators S , obey the Coxeter relations (|3.16[) , (|3,17p and (|3.18p 
for the permutation group S4. These operators depend on the differences of parameters 

Si(u) = Si(ui - u 2 ) ; S 2 (u) = S 2 (u 2 - v x ) ; S 3 (u) = S 3 (ui - v 2 ), 

so that the defining equations can be represented in simpler form 

S fc (-a) Sfc(o) = 1 ; Si(a) S 2 (a + b) Si (6) = S 2 (6) Si (0 + b) S 2 (a) , (5.1) 

Si(a) S 3 (6) = S 3 (6) Si(o) ; S 2 (a) S 3 (a + b) S 2 (6) = S 3 (6) S 2 (a + 6) S 3 (a). (5.2) 
In the theory of quantum integrable systems the cubic relations are known as the star-triangle relations 



There are two evident equalities. The operator S3 (a) differs from the operator Si (a) only by the 
change of variable z\ — > z 2 and therefore these operators commute Si(o) S3 (6) = 83(6) Si (a). 

The equality S 2 (a)S 2 (— a) = 1 is easily checked (for ip = 1) because the operator S 2 reduces to 
the multiplication by a given function. However, if we would have an arbitrary periodic function 
(f(zi, z 2 \ a) in the definition of S 2 (a), then this condition would not fix ip, there are many nontrivial 
periodic functions of zi, z 2 obeying the constraint <f(z±, z 2 ; a)(p(zi, z 2 ; —a) = 1. Only if <p(z\, z 2 ; a) does 
not depend on Zi,z 2 , this condition fixes S 2 up to a constant multiplier <p(a) satisfying the constraint 
ip(a)ip(-a) = 1. 

Let us show that the remaining nontrivial Coxeter relations follow from the elliptic beta integral 
evaluation formula [13]. Remind that the operator S 2 (a) is a multiplication operator 

[S 2 (a) $] (zi,z 2 ) = S (21, 22) • $(zi,z 2 ) ; S a (z 1 , z 2 ) = T(^z 1 ^ z 2 + a + r] + ^) , 

and operator Si (b) is an integral operator 

[Si(6)$](zi,z 2 )= f dzn(z)W b {zi,z)$(z,Z2) ; 
J 

= Y^zJ ' Wb(Z1 ' Z) = 6 T(-2b) T{TZ T Zl ~ b) ■ 

Similarly, 

[S 3 (6)*](^i,^) = f 1 dzfi(z)W b (z 2 ,z)^( Zl ,z). 
Jo 

The first Coxeter relation 

Si (a) S 2 (a + b) Si (6) = S 2 (6) Si (0 + b) S 2 (a) (5.3) 

is equivalent to the following equation for the kernels 
1 

dz n(z) W (£i, z) S a+b (z, z 2 ) W b (z, x) = S b (zi,z 2 ) W a+b (z 1 , x) S a (x, z 2 ) . (5.4) 
The second Coxeter relation 

S 3 (o) S 2 (a + b) S 3 (b) = S 2 (b) S 3 (a + b) S 2 (a) (5.5) 

is equivalent to a similar equation for the kernel 
1 

dzn(z) W a (z 2 ,z) S a+b (z 1 ,z)W b (z,x) = S b (z!, z 2 )W a+b (z 2 ,x) S a (zi,x) . (5.6) 
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This equation can be obtained from the first one after permutation z\ z 2 provided the function 
SaC^i, z 2 ) is symmetric, S a (zi, z 2 ) = S a {z 2 , z%), which is evident in our case. Therefore we have to prove 
the first Coxeter relation only. 

Let us compare the elliptic beta integral [TT] 

J{> y ' v ' k=i i<i<k<e 
where t% H + t§ = 2rj + r, Im(£&) > 0, with relation (15, 6p which has the following explicit form 

dz r^y e r(-2 Q ) r^Tzi-al-rfrz^ + a + i + ^ + S)- 6 ^ r( T x=F^-&) = 



9 17 L 1J 



= T(tzi ^z 2 + b + r] + ^)- T ^_ 2a _ 2b ^ T (Tx TZi-a-b)- T(^x =F z 2 + a + 77 + §) . 
After evident simplifications we arrive at the integral 

/ w o \ ' r ( =FZl z ~ a ) r ("F z 2 =Fz + a + 6 + 7?+§) T(^x ^ z - b) = 
Jo 1 rF^J 

= T(-2a) T(-26) T(2a + 2b + 2?? + r) T(^z 2 T z\ + & + V + §) r(=F«i =F x — a — 6) T(^x ^ z 2 + a + r? + 5) • 
This formula coincides with elliptic beta integral for the following choice of parameters 

ii = z\ — a , t 2 = -z\ - a , t 3 = z 2 + a + b + r] + ^ , 
U = z 2 -a-b + r]+^, t 5 = z 2 - b , t& = -z 2 - b . 

Our normalization ip = 1 (or, more generally, the demand that (p does not depend on zi,z 2 ) is 
not crucial for the obtained result since it will be valid for all (existing) nontrivial periodic functions 
satisfying the constraint 

tpi(z 1} z; a)(p 2 {z, z 2 ;a + b)ipi(z, x; b) = tp 2 (zi,z 2 ; &)</?i(zi,a;; a + b)ip 2 (x, z 2 ; a), 

where <pi(zi, z; a) is associated with the weight W a (z\,z) and (p 2 (zi,z;a) is associated with S a (zi,z). 
Therefore this functional freedom should be fixed in some other way. 
The remaining Coxeter relations 

Si(o)Si(-a) = l; S 3 (a)S 3 (-o) = l, 
are equivalent to the integral relation 

1 1 

dz fi(z)W a (z u z)W- a (z, z 2 ) = [S(z 1 -z 2 ) + 8{ Zl + z 2 )} . (5.7) 

2//(zi) 

After substitution of the explicit expressions for W^-kernels, we obtain the following integral identity 
f r(T^l TZ-a). T(tz 2 T z + a) = T ^a)T{^2 Zl ) { ^ _ ^ + ^ + ^ ] ^ 

which was established in |19j . A more detailed discussion of the intertwining operator properties is given 
in the following two sections. 
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6 Reduction to the finite-dimensional case 

The intertwining operator for integer values of 21 was constructed in [24] . In [23] it was used for building 
solutions of the YBE along the same lines as described above using the realization of Si^a) as difference 
operators of a finite order. In particular, the operator Si (a) had the following form (see the beginning 
of Sect. 5.1 in [23]): 

Q („\ - ^iarjx+ni^ T(2r]x) ^ [2k - X - a] [ j - X - a] [j - a] J a _ 2 k)d x (a Q\ 

r( 2?? (x + a)) ^ \-x-a] jV[ 3 + l- x ][ J + l] e ' W 

where [x] = 6i{2t]x\t). This operator is not well defined unless the infinite series terminates. In 
particular, the action of this operator on a meromorphic function of x yields in general a diverging 
series. Therefore we assume, as in [23], that a = 2t + 1 is an integer. 
Let us denote 

q = e 4ni \ p = e 2n[r , w = e ~ 2nirix , t = e~ 27Tiria . 
Then, using the short theta function, we can write 

ax.ai T(w~ 2 ;p,q) 



Si(a)/H 



TCt 2 w 2 \p,i 



OO 9. 9fc \ k—1 



,8(t 2 w 2 q 2k ;p) ^ 9(t 2 w 2 qi,t 2 qi;p) 
where 9{ti, . . . ,t k ;p) = 6(h;p) . . .6(t k ;p) and 



T{x;p,q) = Y\ 



1 — x 1 p t+l q3+ 1 



1 — xp l q3 
i,j=0 H 

is the elliptic gamma function in the multiplicative notation. 

Introduce now an integral transformation operator acting in the space of z — > z^ 1 invariant functions, 
f(z~ 1 ) = f(z), 

m)f]H = W ^—J -9(z,w,t)^^^lf(z), 
T{t 2 ;p,q) J c z T{z ±2 ;p,q) 

where g(z,w,t) is some fixed function to be determined from the comparison with operator Si (a) (|6.9p 
and 

(p;p)oo(q;q)oc 
K = 'a — '■ • 

47T1 

The contour of integration C is chosen in such a way that it separates geometric sequences of the 
integrand poles converging to zero z = from their reciprocals. In particular, if \tw\, \t/w\ < 1, 
g(z,w,t) has not poles, and f(z) is an analytical function, then one can choose C = T, the unit circle 
with positive orientation. 

Now we take \twq k \ > 1, k = 0,...,N, \twq N+1 \ < 1, \twp\ < 1, \t/w\ < 1, and the contour of 
integration C as a deformation of T including into itself the poles at z = twq k , k = 0, ... ,N, lying 
outside T and excluding the poles at z = t~ 1 w~ 1 q~ k , k = 0, . . . , JV, which enter T. Now we pull C to T 
and pick up the residues at z ±x = twq k , k = 0, . . . , N. This yields the formula 

wiw—f Hm (z-t m k ) T{t zTiT' )P \ q) 

1 KJJiK ' T(t 2 ;p,q) ^ z^twgk^ qj zT(z± 2 ;p,q) 

X (g(twq k ,W,t) +g(t- 1 w- 1 q- k ,w,t)^j f(twq k ) 
k f dz . T{tw ±l z ±l ]p,q) 
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Taking into account the relation 



lim (1 - z)T(z;p,q) 



*-n ' ' {p■,p)oo{q■,q) o , 

we find 

twq k \ . _ 1 . 1 



(twq k \ 
1 I T(tw z~ l ;p, q) 
z J 



z^twq* V z J ' (P;p)oo{q;q)oo0(q k ,...,q 1 ;p)' 

As a result, we obtain 

T(u;- 2 ;p,g) ^ g(twq k ,w,t) + g(t- l w- 1 q~ k ,w,t) 

WWW = r ( t-^ ;Rff )g 2 

x +-4k ~2k ~k 2 8(t 2w2 q 2k 'iP) TT 9(t 2 qi,t 2 w 2 qi;p) k 
Xt w 1 e(t 2 w 2 ;p) f} Q 9(qi+\w 2 q^;p) ntWq } 

dz r(tty ±1 z ±1 ;p, 



Now we take the limit i 2 — > q~ -1 , which means a = 2£ + 1. From the previous conditions 
|iwg' Ar+1 1 < 1 and \t/w\ < 1, we find the constraints \q~ e ~ 1 / 2 \ < |io| < \q e ~ N ~ 1 / 2 \. 

Since 1/T(i 2 ;p, q) = l/r(g~ 2£ ~ 1 ;p, q) = 0, and the integral J* T is finite, we see that the second term 
disappears with the final result 



r(w sr^ 9{wq 2 ,w,q 2 ) + x g + 2 

^)/]H = r( ^ +lu; - 2 . p>(z) g 2 

9 9 %- 2 ^-% 2 ;p) JLA e( g ?+ 1 ,«;V +1 ;p) 

Equating this expression with the discrete series operator Si (a), we come to the relation 

1 2 2 

- (g(z, w, t) + g(z~ l ,w, t)) t~ 4k w~ 2k q~ k2 = q^ + ^~t~ 2k , z = twq k . 

Its solution 

ax, a 2 r(k>gz) 2 - (log to) 21 

L l °sq 

is invariant under the transformation z — > l/z, as wanted. Let us denote z = e 2nir,u . Then we can write 

g ( z = e 2nivu ,w = e~ 27rivx ,t = e" 2 ^) = e ^(« 2 -* 2 ). 
Finally, our integral operator gets a nice form 

m)f]{w)= K f ^-^ yy-y) ^ 

L Wl/U ; r(t2;p, 9 ) 7 C z T(z± 2 ;p,q) M ; 

with the connection between the parameters u, x, a and 2, w, t indicated above. Since the kernel function 
is not analytical in z the contour integral used in the definition of B(t) is formal - to give a proper 
meaning to this operator one has to pass to real integral by a change of variables. 
Let us replace in B(t) 

u x 
u — > , x 



2-K7] ' 2ttt] ' 
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after which we have z = e lu and w = e lx . We change also the integration measure 

* i r* 1 = 2,'' 1 dU 



h z Jo J-i Vl - U 2 

where U = cos u. After these changes our integral operator takes the form 

lBim * ) = ^m¥^Lvrrw e ' /( " ) - 

For t = e~ 27Tlva this operator coincides with the Si(a)-integral operator derived earlier. 
Let us compute the result of acting by the product B{t)B{t~ l ) on a function, 

-\\*\f \ &p) 2 (r,q) 2 f 1 dv Kv 2 -* 2 ) r(tA ±v±x )-. P] q) 



[B(t)B(t- l )f](x) = J™ , —lL^ e ^ 



27r) 2 T(t± 2 ;p,q) J_ x VT^V 2 T(e± 2 ™;p,q) 

1 dU ii^hr(t- l e i( - ±u±x ^;p,q) p . . 

e 47rr > — =T7 Ttv s /(m), 



l-l Vl^TP T(e ±2iu ;p, 
where V = cos v. Changing the order of integrations we obtain 



-j.,,-], , _ (p;p) (<?;<?) I du e 



dV T(te i( - ±v±x "> ; p, g)r(t- 1 e i ( ±u±:r ');p, g) 



4\/rF r(e ±2 ™;p,g) 
Finally, denote X = cos x and use the identity established in [19j[2T] 



(2vr) 2 J_ x T(e± 2iu ;p,q) ^/^^ 



= T(t ±2 , e ±2 ^;p,q)Vl-V 2 5(V - V), (6.11) 
where U = cos u, V = cos 0, V = cos <fi' . As a result, we obtain 

dUe *™> 6(U - X)f(u) = f(x), or B(t)B(t~ x ) = 1, 

as required. 

Comparing 5(f) with the Si(6)-intertwining operator we see that 

[S 1 (6)<I>](x) = [ J B(e- 2 ^)/](e 2 ^), 

where <3?(x) = f(e 2mx ), i.e. these integral operators coincide (the difference is only in the additive or 
multiplicative notations). We conclude that the discrete intertwining operator of [23] is a special limiting 
case of our intertwining operator Si (b). 



7 Connection with the elliptic Fourier transform 

The Bailey chains technique is a well known tool for generating infinite sequences of identities for plain 
and q-hypergeometric series [16] . They are very useful for proving Roger s-Ramanuj an identities needed 
for solving 2d statistical mechanics models [1]. An integral analogue of the Bailey chains was discovered 
in [17] directly at the elliptic level using a new universal integral transform for functions depending on one 
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parameter t. As shown in |28j the inverse of this transform is equivalent to the inversion t — > t~ l which 
resembles the Fourier transform. Let us describe key ingredients of this elliptic Fourier transformation 
(EFT) technique. 

For a given function a(z,t) analytic near the unit circle z € T define the integral transformation 

P(w,t) = M{t) wz a{z,t) := — / — a(z,t) — . (7.1) 

47ri Jj i [t , z ;p,q) z 

In [17] the functions related in this way were said to form an integral elliptic Bailey pair with respect 
to the parameter t (this particular renormalized form of the definition was presented in [28J ) . 

An integral analogue of the Bailey lemma provides a method to produce infinite chains of Bailey 
pairs from a given germ pair. Namely, suppose that a(z, t) and f3(w, t) form an integral elliptic Bailey 
pair with respect to the parameter t. Then the new functions 

a'(w,st) = D(s,u) w a(w,t), D(s,u) w = T(^/pqs~ 1 u ±1 w ±1 ;p,q), (7-2) 
(3'(w, st) = D{t- l ,u) w M(s) wx D(st, u) x f3(x, t) (7.3) 

form an integral elliptic Bailey pair with respect to the parameter st. 
The proof of this statement is easy. Indeed, the demand 

j3'{w,st) = M(st) wz a'(z,st) 

is equivalent to the equality 

D{r x ,u) w M{s) wx D(st,u) x M{t) xz a(z,t) = M(st) wz D{s,u) z a(z,t). 

Since D(i _1 ,u) w D(t,u) w = 1 due to the reflection equation for the elliptic gamma function, we can 
rewrite it as an operator identity 

M(s) wx D(st, u) x M(t) xz = D(t, u) w M(st) wz D(s, u) z . (7.4) 

After substitution of the expressions for M and D-operators one can easily see that it is equivalent to 
the elliptic beta integral evaluation formula which proves the statement. 

Relation (|7.4p is nothing else that the star-triangle relation written in the operator form and it is 
equivalent to the star-triangle relation for the Boltzmann weights D(t,u) w presented in |29| . Equality 
(17. 4p was explicitly presented in [11] in the matrix form as relation (6.5) (it corresponds to the finite 
dimensional reduction of the intertwining operator discussed in the previous section). So, we may 
conclude that the Bailey lemma established in [17] is equivalent to the star-triangle relation for particular 
elliptic hyper geometric Boltzmann weights. 

In [28] it Wcis shown that under certain restrictions the integral transformation ]\/[(t^ wz goes to its 
inverse after the change t — > t~ x . Using real variable integrals as in the previous section, one can write 
M(t -1 )M (t) = 1 |19p27j . where we omit the integral (matrix) operator indices for simplicity. Comparing 
the operator M(t) with the intertwining operator of the previous section we see that 

M{t) wy = e wix2/r >B(t)e- wiu2/ri , (7.5) 

where w = e 2 ™^ and y = e 2mu l' n . So, these two operators differ by a very simple similarity transfor- 
mation. 

For our further purposes we redefine our basic intertwining operator 



Si (a) := M(e- 2 ™U, 
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where we assume that w = e 27T1Z1 and x = e 2mz with variables Z\ and z used for defining Si (a) in 
the previous notation. After performing appropriate similarity transformation for the Sklyanin algebra 
generators, we have 

Si • S J (z, g) = S J (z, -g) ■ Si , j = 0, . . . , 3, 
where g = r](2£ + 1) is an arbitrary parameter, and 

S^z,g) = e^^^ ongmal e-- 2 ^ 
(i)^0 i+ i(r?) 



e •» 



'2z\ 



l (2z-g + r j )- e nd - 9 j+1 {-2z -g + rj)- e~ vd -e~~ z ~ . (7.6) 



Evidently, 82(a) = D(e~ 27Via ,e mZl ) e 2-Kiz 2 - Therefore one easily sees that the operators entering the 
integral Bailey lemma D and M define elementary transposition operators satisfying inversion relations 
5*2 (a) S 2 (—a) = 1, Si(a)Si(—a) = 1 and the basic relation (|7.4|) is equivalent to the Coxeter relation 

Si(a)S 2 (a + 6)Si(6) = S 2 (6)Si(a + 6)S 2 (o). (7.7) 

Defining 83(a) as an EFT transformation operator for variable z 2 we obtain another generator of the 
group 64. So, the algebraic relations for the Bailey lemma ingredients appear to be equivalent to the 
Coxeter relations for permutation operators. This fact was not understood in the original paper |17j . 
Acting by the operator identity (|7.4|) (or Coxeter relation (|7.7|) ) onto a delta- function located at one 
point one obtains the star-triangle relation in the form presented in [29] whose relation to the integral 
Bailey transformation was briefly remarked in 



8 Uniqueness of solutions and the elliptic modular double 

The concept of elliptic modular double introduced in [19] allows us to fix uniquely the form of the 
64-permutation group generators. 

The general solution of Yang-Baxter equation we obtain 

Mi 2 (u) = P12S2CU1 - v 2 ) Si(«i - vi) S 3 (u 2 - v 2 ) S 2 (u 2 - vi) 

is symmetric in the parameters p and q since we assume that U\ 2 ,v\ 2 are independent variables (or, 
equivalently, u and g = r/(2l+ 1) should be considered as independent variables). This operator has the 
following explicit action on the functions of two variables 

[M 12 (u)f]( Zl ,z 2 ) = ^|4|^ r(^g^ 1 4 1 e 2 ^-^);p ;g ) (8.1) 

r^Tri^i-uO^il,^ g&rifa-uaJ^ljfH ^ e 2,ri(«a-«i) a .±l I/ ±l ; p ) q ) dx dy 

r ^ e 4m(«i-ui) } e ^\(v 2 -u 2 ) j x ±2 j y ±2. p ^ g ) ^ X,V '~^~^' 

Because of the symmetry in p and q, we have not one but two RLL-intertwining relations. The 
second one is obtained from (|3.7p simply by permuting 2r/ and r (or p and q): 

R 12 (u - v) Lf uh {u) a 3 L^iv) = L d 2 oub (v) a 3 Lf oub (u) R 12 (u - v) . (8.2) 



Remind that after our similarity transformation the initial L-operator is given by the expression (|2.2p 
where Wj(u) = 9j + \(u + r]\T)/9j + i(r]\T) and Sklyanin generators S J have the form (|7.6p . The operator 
L doub (u) is obtained from L(n) simply by permuting 2r\ and r. This means that L doub (u) also has the 
form (|2.2p . but now Wj(u) = 9j + \(u + t /2\2rj) / 9j + \(T /2\2rj) and the operators S J should be replaced by 



1. 



-9 j+1 (-2z-g + T/2\2ri)-e-^ a -e~^ z \ (8.3) 



2?ri ,2 
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where (/-parameter is the same arbitrary parameter as in (|7,6p . 

Such a direct product of two Sklyanin algebras was introduced in [19] under the name an elliptic 
modular double. It represents an elliptic analogues of Faddeev's modular double [18J. 

The additional symmetry in p and q in our operators S«(a) and the R-operator originates from 
a particular choice of the arbitrary periodic phase factors <p(zi,z 2 ) emerging from the solutions of 
corresponding finite-difference equations. We can invert the logic and demand from the beginning 
existence of the elliptic modular double (|8.2p . Then we can repeat the same considerations as before 
and get the same solutions of the finite-difference equations, but now the phase factors are restricted by 
a different set of periodicity requirements 

(p(zi +t,z 2 ) = tp(z!,z 2 + t) = ip(zi +t/2,z 2 +t/2) = ip(z ll z 2 ). 

Finally we notice that all our operators are meromorphic functions of e 2wlZl and e 2mZ2 . Therefore it 
is natural to demand that we are actually working in the space of functions meromorphic in these two 
exponential variables. This demand imposes additional constraints on the periodic functions we are 
discussing: 

<fi(Z! + 1,Z 2 ) = lf(z 1 ,Z 2 + 1) = ip(zi,Z 2 ). 

Taken altogether, these three conditions yield (p(zi,z 2 ) = const since there do not exist nontrivial triply 
periodic functions. Thus (p may depend only on the parameters tii2jt>i,2> an d we have chosen the 
normalization <p = 1. 

Thus the elliptic modular double (i.e., two RLL-relations for a given R-operator) fixes the intertwin- 
ing operators Sj uniquely (up to the multiplication by a constant). Note that the generators S 3 and S J 
do not commute with each other, however they satisfy some simple anticommutation relations [19] . As 
described also in [19J, there exists a second modular double obtained by employing a modular trans- 
formation of theta functions and the modified elliptic gamma function |14| . but we do not consider it 
here. 

9 Conclusion 

In this paper we have merged two constructions from the theory of quantum integrable systems and the 
theory of special functions. One construction is a specific approach to building YBE solutions developed 
in pOH23] . It is based on the realizations of the permutation group generators by various operators 
acting in the functional spaces and it directly leads to the factorized form of the R-matrix as products 
of elementary transpositions. Another construction is the elliptic beta integral [13] and its various 
consequences formulated as elliptic Fourier transformation and integral Bailey lemma [17j . It formed a 
base for the development of the theory of a principally new and very powerful class of special functions 
- elliptic hyper geometric integrals |lipi4j. As a result of our considerations, both fields have benefited 
and mutually enriched each other. The most general R-matrix at the elliptic level appeared to be defined 
by an integral operator with an elliptic hypergeometric kernel and algebraic properties of the integral 
Bailey lemma ingredients got a natural interpretation as Coxeter relations for the permutation group 
generators. Moreover the key integral operator defining the elliptic Fourier transformation appeared 
to be an intertwining operator for the Sklyanin algebra. The general construction shows that YBE is 
nothing else than a particular word identity in the group algebra for the symmetric group ©4, whose 
generators are realized as integral operators. 

Our results can be applied to all known forms of YBE [7J. In particular, a generalization of our 
construction to root systems is relatively straightforward due to the abundance of elliptic beta integrals 
on root systems |llj and corresponding elliptic Fourier transformations [28j . In the rational case the 
general R-operator for ^4 n -root system was constructed in [21,22j. Star-triangle and star-star type rela- 
tions for the root systems following from the elliptic hypergeometric integral identities were considered 

in [2gE2]!. 
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In |31| , Faddeev and Volkov constructed a solution of YBE at the g-hypergeometric level with the 
help of the pentagon relation for noncompact quantum dilogarithms. A generalization of this model has 
been found in [27] and a question was posed - is it still related to pentagon relation and does there exist 
an elliptic analogue of the latter ? As shown in [23] , the method used in the present paper works at the 
compact g-hypergeometric level. The non-compact situation can be treated as well after appropriate 
change of (/-exponentials by noncompact quantum dilogarithms. However, it is easy to degenerate our 
elliptic results to both compact and non-compact g-levels. From our analysis of the elliptic level case it 
is not clear which relation can be taken as a direct elliptic analogue of the pentagon relation. Instead of 
a potential five-term relation, the key role is played by the hexagon relation (|7.4p emerging in the theory 
of elliptic Fourier transformation [T7] and defining the Coxeter relation for permutation operators (or 
the star-triangle relation) in integrable models of statistical mechanics. 

We would like to stress that the spin variable I in our analysis takes continuous values. Therefore, 
strictly speaking, we deal not with the discrete Ising-type models, but with two-dimensional quantum 
field theories. In this context the Yang-Baxter equation can be interpreted as a condition of factorizing 
the iV-body 5-matrix to the product of two-body scattering matrices [7J. 

Let us discuss briefly application of our results to four dimensional super symmetric gauge field 
theories. The key discovery of [32] consists in the fact that superconformal indices of these theories 
are described by the elliptic hyper geometric integrals. For example, the elliptic beta integral gets a 
remarkable interpretation as the direct indication on the confinement phenomenon in the simplest four- 
dimensional supersymmetric quantum chromodynamics. 

A relation between four-dimensional Nekrasov instanton partition function and conformal blocks in 
two-dimensional Liouville field theory was empirically discovered in [33] ■ In [31], superconformal indices 
of Ad N = 2 supersymmetric field theories were tied to correlation functions of 2d topological field 
theories. Our results are relevant to a different type of Ad/ 2d correspondence discovered in [27], where 
Ad superconformal indices coincide with partition functions of integrable models of 2d spin systems. 

Seiberg duality is a special electric-magnetic duality of four dimensional supersymmetric non-abelian 
gauge field theories. In the language of elliptic hypergeometric integrals there are two qualitatively 
different situations. When the dual theory confines corresponding superconformal index is identical to 
some elliptic beta integral on a root system or, from the statistical mechanics point of view, to star- 
triangle relation for multicomponent spin systems. When the dual theory is a nontrivial interacting 
field theory, one deals with symmetry transformations for integrals equivalent to the star-star relations 
in statistical mechanics [271. 



Since superconformal indices coincide with the statistical sums of elementary cells for plane spin sys- 
tems, the Seiberg duality transformations represent the Kramers- Wannier type duality transformations 
for corresponding new 2d integrable spin systems. As indicated in |35j sequential integral transforma- 
tions following from the Bailey lemma define superconformal indices of quiver gauge theories. This 
procedure corresponds to building full two-dimensional lattice spin systems as prescribed in the theory 
of quantum integrable systems. 
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10 Appendix 

In this Appendix we collect some useful formulae. The standard infinite g-product is defined as 



+00 



(x;q)oo = ~q k -x) ; qeC, \q\<l. 



k=0 

The general theta-function with characteristics has the form 



nGZ 

We use four standard theta-functions 



where p = e 27nr and 



n£Z 

= ip 1 / 8 e~™ ( P ;p)ooe(e 2 ™;p), 

0(t;p) = (t;p)oo(pi _1 ;p)oo, 
9 2 (z\t) = 1>0 {z\t) = e wi ("+l) 2 - . e M n H> 

neZ 

3 (z\t) = e 0fi (z\r) = J] e™ V • e 2 ™ 2 

neZ 

9 4 {z\t) = 0o,i (z|t) = ^e™ V •e 2 " i ("+l)( z +^). 

neZ 

The following identities are used to factorize the Lax matrix and to derive defining equations for the 
operator S2: 

2 9 1 (x + y)9 1 (x-y) = 9 4 (x) 9 3 (y)-9 4 (y)9 3 (x) ; 2 9 4 (x + y) 6 4 {x - y) = h{x) 6 3 (y) + h{v) h{x)\ 

2 6 2 (x + y) 6 2 (x - y ) = 9 3 (x) 9 3 (y ) - h (y) 9 4 (x) ; 2 9 3 (x + y) 9 3 (x - y ) = 9 3 (x) h (y ) + 4 (y) h (x) , 

where 6 3 {z) = 9 3 (z|§) and 4 (z) = 4 (z\^). 

For Im(r) > 0, lm.(rj) > the elliptic gamma function is defined by the double infinite product 

00 1 _ e 27ri(r(n+l)+2r?(m+l)-2) 

n»\T,2n)= II i _ e27ri( ,n +2 , m+2) ■ ( 1(U ) 

n,m=0 

It satisfies equations 

T{z + t\t, 2ri) = 9(e 2wiz ; e 4wir ?) • T(z\t, 2rj) (10.2) 
T(z + 2 V \t, 2rj) = 0{e 2 ™; e 2 ™) ■ T{z\t, 2 V ), (10.3) 

One can evidently write 

9(e 2niz ; e 27riT ) = R(r) • e^i^r), 9(e 2wiz ; e 4 ^) = R(2r?) • e^O^zftri), 
where the constant R(r) does not depend on z: R(r) = — ie 4" • (e 27rlr ; e 2,riT ) 1 . 
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